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I [NTRODUCTION  AND  SUMMARY 


This  semiannual  technical  report  describes  major  aspects  of  the 
technical  support  provided  by  SRI  International  to  the  Electronic  Sys 
terns  Division  (ESD) , Air  Force  Systems  Command  (AFSC) , in  the  develoo- 
ment  of  the  CONUS  414L  Experimental  Radar  System  (ERS) . The  work  was 
performed  during  the  period  1 January  through  30  June  1977  under  Con- 
tract  N00014-75-C-0930.  To  resolve  questions  pertinent  to  risk  areas 
in  the  ERS,  SRI  has  used  the  Wide  Aperture  Research  Facility  (WARF) , an 
existing  over-the-horizon  backscatter  (OTH-B)  HF  radar.  WARF  is  an  SRI- 
operated,  ESD/ONR-sponsored  facility  capable  of  simulating  ESD's  ERS 
functions  in  almost  all  respects,  although  on  a smaller  scale,  WARF 
serves  as  a test-bed  for  newly  developed  OTH-B  techniques  related  to  the 
ERS  procurement. 

The  CONUS  414L  Experimental  Radar  System  currently  being  procured 
is  an  over-the-horizon  HF  backscatter  radar  designed  to  detect  and  track 
all  aircraft  passing  through  the  North  Atlantic  air  corridors.  Require- 
ments for  extended  range  coverage  and  broad  azimuthal  sectors  dictate 
the  use  of  a large  energy  product,  and  time-sequential  range-changing  and 
azimuth-scanning  capability.  Normal  commercial  traffic  passing  through 
the  radar  coverage  requires  that  a large  number  of  tracks  be  maintainable 
at  all  times.  Facilities  for  real-time  correlation  and  identification 
will  be  implemented  at  the  radar  receive  site.  The  azimuth  sidelobe 
levels  of  the  receive  antenna  system  are  presently  specified  at  -50  dE 
as  a means  for  assuring  ERS  operation  near  the  auroral  zone  without 


’'Formerly  S ford  Research  Institute. 

^Contract  N00014-75-C-0930  was  issued  by  the  Office  of  Naval  Research  as 
a continuation  of  Contract  N00014-70-C-0413.  Its  sponsors  (each  for 
different  aspects  of  OTH  technology)  are  ESD  (AFSC),  ONR,  and  ARPA. 

The  Scientific  Officer  for  the  overall  contract  is  Mr.  John  J . Kane, 

ONR  [Code  461(FP)],  telephone  (202)  692-4203. 
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reception  of  intolerable  levels  of  auroral  clutter.  The  ERS  system 
performance  test  is  currently  scheduled  to  begin  in  FY80.  Conclusions 
drawn  from  the  ERS  design  will  support  the  transition  to  the  fully 
operational  radar  system  (RS) . 

During  this  reporting  period,  work  on  the  project  at  SRI  has  been 
directed  toward  the  following: 

• Use  of  a subclutter  visibility  (SCV)  measure-.ient  technique  for 
collecting  round-the-clock  data  on  potential  radar  sensitivity 
for  any  available  choice  of  radar  parameters. 

• Development  of  a technique  for  the  suppression  of  radio  fre- 
quency interference  (RFI) . 

• Implementation  of  a model  for  OTH-B  radar  data  to  test  the 
sensitivity  of  the  ERS  display  of  Doppler  (nested  range)  vs. 
time  and  its  detection/tracking  (D/T)  systems. 

• Development  of  a generalized  sidelobe  canceller  (CSC)  beam- 
forming  processor. 

• Stu.  ' f spread-clutter  effects  on  adaptive  beamformers  and  pre- 
liminary testing  of  the  GSC  against  azimuthally  spread  clutter 
generated  by  mis-steering  the  WARE  transmit  beam. 

• Investigation  of  architectures  for  hardware  implementation  of  an 
adaptive  bearaforming  processor  like  the  GSC. 

• Provision  of  direct  support  to  the  ESD  414L  SPO  through  special 
on-request  studies,  visits,  and  consultations  concerning  the  ERS. 

The  status  of  each  effort  is  briefly  summarized  below. 

The  round-the-clock  measurements  involve  a year-long  data  collection 
effort  now  near  its  midpoint.  Analysis  of  the  soundings  has  been  made 
with  the  recognition  that  the  current  operating  concept  for  the  ERS  will 
provide  a single  radar  fence  (barrier)  that  moves  in  range  as  propagation 
conditions  change.  First  analyses  indicate  that: 

• Azimuthal  variations  of  SCV  are  small  over  most  of  the  diurnal 
cycle,  so  the  best  operating  parameters  for  an  azimuth-scanning 
radar  may  be  obtained  by  sampling  only  a few  azimuth-  1 bearings. 
The  major  exception  to  this  result  will  occur  at  ERS  during 
certain  periods  of  auroral  activity  when  it  is  expected  that  one 
edge  of  the  azimuthal  coverage  sector  will  be  up  against  the  edge 
of  intense  auroral  scatter. 

• Usable  barrier  depth  is  limited  to  ranges  less  than  500  nmi.  Two 
24-hour  tests  show  that  the  outer  half  of  the  500-nmi  barrier  is 
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substantially  less  sensitive  (often  as  much  as  20  dB)  than  the 
inner  half.  This  state  of  affairs  should  generally  hold  for 
the  ERS  environment  as  well. 

• Predicted  nighttime  radar  sensitivities  are  higher  than  those 
measured  at  WARF  during  24-hour  aircraft  surveillance  tests.  A 
number  of  explanations  are  possible,  one  of  which  is  that  there 
has  been  a reduction  at  the  radar  of  atmospheric  and  man-made 
noise  originating  from  the  east.  This  reduction  is  most  likely 
the  result  cf  installation  of  256  twin-whip  endfire  receiving 
pairs  (WERP)  in  place  of  the  256  single  vertical  whips  used 
during  the  24-hour  aircraft  surveillance  tests. 

• Large  and  frequent  inward  or  outward  movements  of  the  barrier 
will  be  necessary  if  higher  radar  sensitivity  is  to  be  maintained. 
A tradeoff  between  sensitivity  and  barrier  movements  will  be  re- 
quired to  avoid  missing  targets  through  discontinuities  in  bar- 
rier coverage. 

• Existing  off-line  algorithms  designed  as  automatic  aids  for  se- 
lecting the  best  radar  operating  parameters  perform  consistently 
well,  and  could  become  valuable  aids  in  on-line  propagation 
management,  even  though  they  do  not  currently  account  for  multi- 
path,  auroral  clutter,  or  frequent  large  barrier  movements. 

RFI  is  an  anticipated  problem  for  the  ERS.  MITRE  expects  a typical 
RFI  environment  of  several  (3  to  5)  carrier  sources,  each  having  approxi- 
mately a 20-dB  SNR  and  leading  to  SNR  degradations  of  20  dB  for  several 
target  speeds.  An  SRi-developed  RFI-suppression  algorithm  known  as 
BRISA  can  be  applied  directly  to  scalar  amplitude  data  after  coherent 
processing  to  reduce  or  eliminate  almost  all  RFI  signal  energy.  Basic- 
ally, BRISA  seeks  by  efficient  computer  processing  to  eliminate  the  bias 
introduced  by  the  constant  level  of  RFI  containing  carrier  components. 

The  expected  20-dB  SNR  degradations  for  several  target  speeds  should  be 
reducible  to  about  3 dB  for  each  of  these  speeds  through  application  of 
BRISA.  Coupling  of  BRISA  to  the  frequency-stepping  technique  proposed 
by  MITRE  should  significantly  improve  ERS  performance,  especially  in  the 
typical  nighttime  environment.  In  this  report  we  demonstrate  the  present 
capabilities  of  BRISA  in  processing  WARF  data. 


BRISA:  Bias-Removal  Interference-Suppression  Algorithm, 
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An  SRI-developed  analytical  model  that  simulates  radar  target  and 
noiP.e  data  provides  a controlled  environment  for  the  test,  evaluation, 
and  comparison  of  the  WARF  animated  display  system  and  the  ERS  grayscale 
display  format.  The  modnl  provides  for  awell-to-dwell  SNR  fluctuations, 
target  signature  changes  Including  range  and/or  Doppler  spreading  of 
target  signatures,  and  randomly  occurring  meteors  at  a prespecified  ex- 
pected rate  of  occurrence.  Not  included  are  long-term  fading  over  many 
dwells  and  signature  variation  within  an  established  track. 

An  initial  test  of  the  analytical  model  was  directed  toward  deter- 
mining the  detection  sensitivity  of  the  ERS  display.  Earlier  tests  done 
elsewhere  evaluated  display  sensitivity  with  nonfluctuating  CW  signals 
in  noise  and  used  progressively  stepped  SNR  levels.  The  model  tests  at 
SRI  sought  to  reduce  the  element  of  bia?  introduced  by  the  monotonic 
steps  of  the  earlier  tests.  A display  of  range  (nested  Doppler)  versus 
time  containing  targets  of  0 to  20  dB  SNR  was  generated.  Our  analysis 
of  the  tracks  within  the  simulated  ERS  display  data  yields  an  estimated 
display  sensitivity  of  between  6 and  9 dB  as  the  SNR  required  to  obtain 
a 507«  detection  rate.  This  figure  is  representative  of  manual  detection 
results;  use  of  the  SEDAT*  algorithm  should  lower  this  threshold  somewhat. 
In  any  event,  the  manual  detection  threshold  as  determined  by  this  model 
of  the  ERS  display  is  several  dB  higher  than  that  determined  by  the  tests 
involving  nonfluctuating  CW  signals. 

The  University  of  Colorado  under  subcontract  to  SRI  has  developed  a 
new  adaptive  noise-cancelling  scheme  referred  to  as  a Generalized  Sidelobe 
Canceller  (GSC) . A most  promising  aspect  of  this  processor  structure  is 
that  it  permits  the  realization  of  a fairly  general  class  of  adap.ive 
beamforming  algorithms  with  a single  unifying  analytical  approach.  A 
version  of  this  processor  has  been  implemented  in  software  at  SRI  for 
off-line  processing  of  WARF  data  and  for  preliminary  performance  evalua- 
tion. The  GSC  is  important  because  with  further  development  it  could 


*SEDAT:  Sequential  Detection  and  Tracking. 
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provide  a fallback  method  for  rejection  of  auroral  clutte.  at  the  ERS 
site  or  for  tha  eventual  R3. 

Application  of  adaptive  beamforming  techniques  to  HF  radar  is  at- 
tractive because  of  their  potential  for  rejecting  spread-azimuth  and 
spread-Doppler  clutter  such  as  that  produced  by  scattering  from  luroral 
irregularities.  A numb< r of  alternatives  have  been  considered  for  using 
the  WARF  adaptive  beamforming  capabilities  to  test  methods  of  rejecting 
spread-azimuth  interference  on  an  experimental  basis  (in  contrast  to 
computer  simulation  studies  that  cannot  reasonably  account  for  various 
aspects  of  the  HF  environment  and  the  constraints  on  HF  equipment) . Op- 
tions include  generation  of  azimuthal  and  Doppler  spreading  by  (1)  using 
a separate  transmitter  to  transmit  a 2-hop  signal  that  is  received  as 
clutter  in  the  sidelobes  of  the  receive  array,  (2)  mis-steering  the  WARF 
transmit  beam  to  generate  a clutter  distribution  with  significant  out--f- 
beam  components,  (3)  collection  of  scatter  from  spread-F  ionization,  or 
(4)  use  of  Ionospheric  heating  to  generate  field-aligned  irregularities 
similar  to  those  found  in  auroral  regions.  Anticipated  frequency  manage- 
ment difficulties  have  caused  us  to  reject  the  separate  transmitter  ap- 
proach. The  mis-steering  option  has  been  trie-’  and  is  discussed  in  this 
reoort.  Observations  of  spread  clutter  due  tc  onospheric  heater-induced 
irregularities  may  be  carried  out  in  FY78. 

In  a preliminary  test  of  the  capability  of  the  GSC  beamformer  to 
reject  spread  clutter,  the  WARF  transmit  beam  was  steered  away  from  the 
desired  target  (the  Los  Lunas,  New  Mexico  repeater)  so  as  to  introduce 
substantial  ground  clutter  in  the  sidelobe  angles  cf  the  receive  beam 
(Option  2 in  the  previous  paragraph).  In  brief,  the  average  signal- to- 
clutter  ratio  (SCR)  improvement  with  the  adaptive  beamformer,  compared 
to  a conventional  beamformer,  operating  in  a clutter  environment  was 
measured  at  8 dB  with  beam  alignment  and  11  dB  with  beam  misalignment. 

The  tests  clearly  demonstrate  t.  at  appreciable  cancellation  of  clutter 

spread  in  azimuth  is  achievable  by  an  adaptive  technique. 

As  a part  of  the  ongoing  analysis  of  adaptive  beamformers  at  SRI,  a 

study  has  been  initiated  to  define  how  an  on-line  beamformer  might 
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eventually  be  best  implemented.  The  GSC  structure  was  used  as  a candi- 
date system  to  obtain  a measure  of  cost,  speed,  and  performance  of  such 
a device  when  integrated  into  the  operation  of  an  OTH  radar.  The  GSC 
structure  led  to  a design  for  a multiplier-oriented  processor  which,  con- 
sidering the  con^jlexity  of  operations  and  the  desire  to  modify  the  GSC 
algorithm  in  real  time,  suggested  a microprocessor-based  control’  if 
the  des.^n  is  ultimately  executed  in  hardware  devices.  The  MOS  cit js  of 
microprocessors  would  be  too  slow,  having  instruction  times  of  at  least 
1 (iS,  so  a bipolar  microprocessor  slice  or  a microprogram  sequencer  might 
be  viable  alternatives.  System  architecture  and  some  aspects  of  hardware 
realization  are  discussed  in  Appendices  A and  B. 

Finally,  with  respect  to  the  provision  of  direct  support  to  the  BSD 
414L  SPO,  SRI  participated  in  conversion  of  the  PRS  concept  to  the  ERS; 
originated  approximately  200  hours  of  sounding  transmissions  for  recep- 
tion by  RADC;  developed  preliminary  design  and  cost  estimates  for  trans- 
ponders usable  in  the  ERS  coverage;  and  developed  a technical  concept 
for  augmenting  the  E7.S  with  an  independently  operated  environmental 
assessment  system,  a minicomputer-based  system  designed  to  improve 
ionospheric  propagation  management  capabilities  at  the  ERS. 


II  STATUS  OF  PRINCIPAL  TECHNICAL  EFFORTS 


A.  Subclutter  Visibility  Measurements 

SRI  has  developed  a method  of  measuring  SCV  from  the  widesweep  back- 
scatter  sounder  (WSBS)  used  for  propagation  management.  This  method  has 
the  advantage  that  radar  sensitivity  can  be  inferred  for  all  possible 
range  gates  and  frequencies,  not  just  those  currently  used  by  the  radar. 
On-line  use  of  this  measurement  technique  during  a WARF  24-hour  aircraft- 
tracking  test  in  August  1975  demonstrated  its  ability  for  optimizing 
the  selection  of  radar  operating  parameters. 

Because  this  WSBS-SCV  technique  is  a direct  measure  of  potential 
radar  sensitivity  for  all  ranges  and  frequencies  it  is  a useful  tool  for 
cnaracterizing  potential  radar  performance.  Such  information  could  be  ^ 
valuable  for  determining  ERS  performance  and  developing  procedures  for 
propagation  management.  For  these  reasons  SRI  has  undertaken  a one-year 
experimental  program  designed  to  collect  24  hours  of  WSBS-SCV  data  ap- 
proximately once  each  month.  In  this  program  a sounding  is  taken  every 
5 minutes.  Three  azimuth  directions  spaced  15°  apart  are  sounded  sequen- 
tially, achieving  a revisit  at  each  azimuth  every  15  minutes.  Each  sound- 
ing covers  from  6 to  27  MHz  in  frequency  and  from  0 to  24  ms  in  time  delay. 

Existing  WARF  computation  facilities  are  inadequate  to  process  and 
display  these  results  in  real  time.  A data-collection  program  was  devised 
that  collects  the  raw  data  and  does  sufficient  real-time  processing  to 
reduce  the  amount  of  data  to  within  manageable  limits.  The  resulting  data 
are  recorded  on  digital  magnetic  tape  for  later  processing  at  SRI  Menlo 
Park. 

During  this  reporting  period,  soundings  from  five  continuous  24-hour 
periods  were  recorded.  Initial  analysis  was  directed  toward  assessing 


References  are  listed  at  the  end  of  this  report. 


7 


the  potential  perforirance  of  an  OTH  radar  operating  in  a three-range-gate, 
three- frequency  format  such  as  that  originally  proposed  for  the  Prototype 
Radar  System  (PRS),  For  each  of  these  range  gates  the  computer  analysis 
involved  a search  to  locate  the  frequency  that  maxim.ized  average  SCV. 

These  frequencies  and  the  resulting  SCVs  were  then  plotted  versus  time. 

In  accordance  with  a change  in  the  proposed  radar  system  design,  this 
analysis  was  later  repeated  with  the  restriction  that  only  the  A,  C,  and 
E high-frequency  radio  bands  were  available  for  radar  use. 

In  mid-February  1977  a design  revision  of  the  414L  program  changed 
the  PRS  to  an  Experimental  Radar  System  (ERS).  The  new  radar  operating 
concept  consists  of  a single  radar  fence  (barrier)  that  moves  in  range 
as  propagation  conditions  change.  SRI ' s SCV  soundings  were  reanalyzed 
to  reflect  this  change  in  operation  concept.  Each  sounding  was  searched 
in  frequency  and  time-delay  to  maximize  the  average  SCV  over  the  range 
extent  of  the  barrier.  Barrier  depths  of  500  and  250  nmi  wete  chosen 
for  rhis  analysis.  Figures  1,  2,  and  3,  respectively,  plot  SCV,  time 
delay  (slant  range),  and  frequency  versus  time  for  the  24-hour  period 
starting  on  12  April  1977.  A 250-nmi  barrier  was  selected  in  the  analysis 
used  to  make  these  plots. 

Thus  far,  analysis  of  the  soundings  has  produced  the  following  find- 
ings: 

(1)  Azimuthal  variations  are  small.  Based  on  analysis  of  two 
24-hour  tests,  it  appears  that  differences  in  operating 
parameters  selected  independently  for  the  three  azimuths 
monitored  were  small  over  most  of  the  diurnal  cycle. 

Thus  it  may  not  be  difficult  to  determine  the  best  operat- 
ing parameters  for  an  azimuth-scanning  radar  by  sampling 
only  a few  azimuth  bearings. 

(2)  Usable  barrier  depth  is  limited  to  ranges  less  than  500  nmi. 

It  appears  that  radar  sensitivity  for  a single  radio  fre- 
quency decreases  rapidly  with  range,  once  past  the  skip 
distance.  Based  on  two  24-hour  tests,  a comparison  of 
average  SCV  for  the  250-  and  500-nmi  barriers  shows  that 
the  outer  half  of  the  500-nmi  barrier  is  substantially 

less  sensitive  (often  as  much  as  20  dE)  chan  the  inner 
half.  For  this  reason  we  have  chosen  to  use  the  250-nmi 
barrier  as  the  basis  for  further  analysis. 
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TIME  — UT 


FIGURE  1 SUBCLUTTER  VISIBILITY  FOR  250-nmi  BARRIER 


(3)  Nighttime  sensitivity  appears  better  than  expected.  As 
seen  in  Figure  2,  the  predicted  nighttime  SCVs  are  quite 
high,  significantly  higher  than  those  measured  from  the 
radar  directly  during  24-hour  aircraft  surveillance  tests 
at  WARF,  It  is  not  yet  known  whether  this  effect  is  real 
or  due  to  complexities  in  the  measurement  of  SCV.  Possible 
explanations  Include  computational  underflow;  algorithm 
difficulties  brought  about  by  the  nighttime  RFI  environ- 
ment; the  presence  of  spread-F  (which  could  reduce  radar 
sensitivity  without  affecting  sounder-inferred  sensitivity); 
or,  finally,  the  use  in  the  WARF  receiving  array  of  TWERPS* 
instead  of  single  array  elements  (which  reduces  atmospheric 
and  man-made  noise  originating  from  the  east).  A 24-hour 
aircraft  surveillance  test  planned  for  FY78  will  yield 


•k 

TWERP:  Twin-whip  endfire  receiving  pair.  Use  of  256  TWERP  elements  in 

the  WARF  receiving  array  yields  a switchable  cardiod  pattern  designed  to 
enhance  reception  from  the  east  or  the  west,  as  selected  by  the  site 
operators. 
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TIME  DELAY 


TIME  — UT — - 

FIGURE  2 VARIATION  IN  TIME  DELAY  (slant  range)  TO  NEARER  BOUNDARY  OF  250-nmi 
BARRIER 


additional  information  on  the  contribution  of  TWERPS  to 
higher  SCV  levels, 

(4)  Strict  maximization  of  SCV  causes  unwanted  gaps  in  radar 
coverage.  The  analysis  programs  choose  an  operating  fre- 
quency and  range  gate  that  maximizes  SCV,  As  seen  in 
Figure  2,  this  criterion  can  result  in  large  fluctuations 
of  the  barrier  position.  To  avoid  missing  targets  through 
discontinuities  in  barrier  coverage,  it  is  desirable  to 
maintain  a fixed  barrier  position  or  one  that  changes  only 
in  small  increments.  Thus  it  appears  that  a tradeoff  is 
often  necessary  between  sensitivity  and  amount  of  barrier 
movement . 

(5)  Automatic  aids  for  propagation  management  appear  achiev- 
able. The  programs  developed  for  off-line  analysis  of 
SCV  soundings  use  automatic  methods  for  selecting  oper- 
ating parameters.  Currently  these  algorithms  do  not 
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FIGURE  3 OPTIMUM  FREQUENCY  FOR  250-nmi  BARRIER 


consider  multipath,  auroral  clutter,  or  frequent  large 
barrier  movements.  Nonetheless,  the  algorithms  do  per- 
form consistently  well  in  selecting  the  best  operating 
parameters,  as  determined  by  post-operating  analysis. 

With  further  modification,  the  algorithms  could  become 
valuable  aids  in  on-line  propagation  management. 

Continuation  of  the  collection  of  round-the-clock  SCV  measurments 
is  planned  for  the  next  six  months.  A technical  report  describing  the 
sounding  system,  methods  of  analysis,  and  early  results  will  be  issued. 


B.  Radio  Frequency  Interference 

To  maximize  sensitivity  of  an  OTH  radar  it  is  essential  to  avoid  or 
suppress  interfering  signals  caused  by  other  radio  transmissions.  In 
normal  operation  of  an  OTH  radar  all  possible  measures  should  be  taken 
to  select  frequencies  that  are  free  of  interference.  Unfortunately,  at 
times  (particularly  winter  nights)  it  will  be  impossible  to  avoid 
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interference  entirely,  and  thus  RFI  suppression  techniques  will  be  nec- 
essary to  maintain  radar  sensitivity.  Also,  since  interference  often 
occurs  unexpectedly,  use  of  a RFI  suppression  algorithm  is  necessary  to 
maintain  operations  until  a new  and  better  operating  frequency  can  be 
selected. 

The  previous  semiannual  report  described  work  conducted  at  SRI  in 
the  investigation  of  spectrum-monitoring  techniques  for  RFI  avoidance, 
and  in  the  evaluation  o.  various  proposed  RFI-suppression  techniques. 

During  the  period  dealt  with  in  the  previous  report,  a promising  new  ap- 
proach to  RFI  suppression  was  conceived.  It  has  now  been  developed  at 
SRI,  and  is  called  the  Bias-Removal  Interference-Suppression  Algorithm 
(BRISA).  A description  of  BRISA  is  presented  here,  along  with  some  ex- 
amples demonstrating  its  ability  to  remove  the  deleterious  effects  of 
RFI.  Perhaps  the  most  remarkable  feature  of  BRISA  is  the  simplicity  with 
which  it  handles  a complex  signal  environment  with  only  a slight  increase 
in  signal-processing  requirements. 

I 

1.  The  Bias-Removal  Interference-Suppression 

Algorithm  (BRISA) 

The  most  troublesome  signal  component  contained  in  most  inter- 
fering radio  transmissions  is  the  coherent  "carrier."  At  the  output  of 
the  radar  processor  this  carrier  is  compressed  in  Doppler,  thus  creating 
a relatively  large  interference  level  that  can  mask  targets  at  the  same 
Doppler.  Cost-effective  methods  have  been  developed  to  spread  this  type 
of  RFI  throughout  all  Doppler  cells  and  thus  minimize  the  average  system 
degradation.  A preferable  approach  would  be  to  reduce  or  eliminate  the 
RFI  signal  energy  itself.  Methods  to  do  this  include  adaptive  beamforming 
and  compressive  filtering.  However,  these  methods  are  costly  in  terms  of 
radar  processor  loading.  The  algorithm  presented  here  is  a relatively 
simple  computation  that  can  eliminate  much  of  the  unwanted  RFI  energy. 

When  used  in  conjunction  with  the  RFI  frequency-stepping  techniques,  this 
new  algorithm  should  be  able  to  eliminate  most  deleterious  effects  of 
RFI. 
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BRISA  is  applied  directly  to  the  scalar  amplitude  data  after 
coherent  processing*  Attributes  of  BRXSA  include  the  following. 

• Requires  only  simple  computations  to  be  performed  on 
radar  data  after  range/Doppler  processing. 

• Affects  only  those  Dopplers  having  range-correlated  RFI. 

• Does  not  affect  range  sidelobes. 

• Can  suppress  RFI  enet^y  by  30  dB  or  more. 

• Can  suppress  multiple  carriers  simultaneously. 

• Can  be  used  in  conjunction  with  phase-  or  frequency- 
stepping RFI-spreadlng  techniques. 

2.  A Short  Description  of  BRISA 

The  effect  of  BRISA  and  the  way  in  which  it  suppresses  inter- 
ference can  best  be  described  with  the  use  of  an  example.  The  plots  of 
Figure  4 show  the  output  from  a single  OTH  radar  dwell.  These  data  were 
collected  at  WARF.  A target  of  opportunity  and  a controlled  RFI  source 


-30  30 


-30  30 


(a)  WARF 

NORMALIZATION 
(RFI/N  - 27  dB) 


DOPPLER  FREQUENCY — Hz 
lb)  ERS 

NORMALIZATION 
(SNR  - 3 dB) 


(c)  BRISA 

(SNR  > 16  dB) 


FIGURE  4 AN  EXAMPLE  OF  INTERFERENCE  SUPPRESSION  USING  BRISA 
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are  contained  within  these  data.  In  each  plot,  amplitude  is  plotted 
versus  Doppler  for  each  of  12  range  resolution  cells.  The  three  plots 
of  Figure  4 each  show  a different  method  of  data  normalization  for  the 
same  radar  dwell. 

Figure  4(a'>  represents  the  normalization  used  for  the  WARF  radar 
devection  display.  One  es.imate  of  noise  is  made  for  the  radar  dwell  and 
all  amplitudes  are  normalized  (divided)  by  this  value.  Very  large  signals 
such  as  the  clutter  near  zero  Doppler  and  the  RFI  at  -27  Hz  are  clipped 
20  dB  above  the  computed  average  noise  value.  This  normalization  has  the 
disadvantage  of  presenting  a nonuniform  background,  but  this  is  partially 
compensated  for  by  the  ability  to  easily  detect  the  presence  of  RFI  energy. 

Figure  4(b)  represents  the  normalization  planned  for  use  with 
ERS.  For  each  Doppler  (column)  an  estimate  of  noise  is  calculated  and 
amplitudes  within  each  column  are  then  normalized  (divided)  by  the  noise 
value  at  that  Doppler.  This  normalization  has  the  advantage  of  whitening 
the  entire  display  field,  but  at  the  cost  of  suppressing  all  information 
on  the  presence  of  RFI.  Wliile  the  presence  of  strong  signals  such  as 
clutter  and  RFI  have  been  eliminated  from  the  display,  this  normalization 
of  course  does  not  improve  target  detectability  at  the  clutter  and  RFI 
Dopplers. 

In  Figure  4(c)  BRISA  has  been  applied  to  the  data  with  the  re- 
sult that  a target  echo  of  approximately  16  dB  signal-to-noise  ratio  (SNR) 
is  now  visible  at  the  RFI  Doppler  frequency.  To  achieve  this  target  de- 
tection, BRISA  has  taken  advantage  of  the  fact  that  while  the  RFI  carrier 
adds  extra  power  at  a particular  Doppler,  it  adds  the  same  arount  of  power 
to  all  range  bins  at  that  Doppler.  Thus,  while  the  average  amplitude  of 
each  range  bin  can  be  greatly  increased  by  an  RFI  signal,  the  standard 
deviation  (a)  of  these  range  bins  is  only  slightly  altered.  Since  detec- 
tion  is  really  limited  by  noise  fluctuation,  (proportional  to  a)  and  not 

*For  Rayleigh-distributed  amplitudes  (a  good  approximation  of  the  noise 
and,  to  a lesser  extent,  clutter  amplitudes)  the  mean  (m)  is  linearly  pro- 
portional to  the  standaru  deviation  (a)  , and  division  by  the  mean  is  there 
fore  appropriate.  For  RFI  (which  is  not  Rayleigh  distributed)  the  mean 
can  be  very  large  while  the  standard  deviation  remains  small.  Hence, 
division  by  the  mean  is  not  the  best  method  of  optimizing  SNR. 
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by  Its  average  value,  a target  is  detectable  despite  the  presence  of 
this  RFI. 

BRISA  seeks  to  eliminate,  before  normalization,  the  bias  intro- 
duced by  the  constant  level  of  RFI.  The  way  the  bias  is  removed  can  be 
demonstrated  by  the  example  of  Figure  5.  Here  the  radar  dwell  of  Fig- 
ure 5 has  been  reprocessed.  For  each  Doppler  resolution  cell  in  Figure 
5(a)  the  rms  amplitude  of  the  12  range  cells  is  plotted  on  a logarithmic 
scale.  Also  plotted  is  the  peak  amplitude  of  the  12  range  cells.  Noise, 
clutter,  and  RFI  are  distributed  in  range,  so  the  rms  plot  is  appropriate 
to  estimate  their  strength.  On  the  other  hand,  a target  is  confined  in 
range  and  here  the  peak  plot  is  appropriate  for  measuring  its  strength. 
Hence,  both  rms  Jind  peak  plots  are  needed  to  determine  such  quantities  as 
SNR  and  slgnal-to-lnterfereuce  ratio.  In  Figure  5(b)  the  ratio  of  the 
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FIGURE  5 DOPPLER  SPECTRA  AND  SIGNAL  VARIATION  FOR  A RADAR  DWELL 
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mean  (m)  to  the  standard  deviation  (a)  of  the  amplitudes  for  the  twelve 
ranges  in  each  Doppler  resolution  cell  has  been  plotted.  For  Rayleigh- 
distributed  values  (such  as  noise  and,  to  a lesser  degree,  clutter)  the 
ratio  of  m/a  should  be  about  2.  It  can  be  seen  from  Figure  5(b)  that 
this  ratio  holds  approximately  true  for  Dopplers  having  both  clutter  and 
noise.  For  Dopplers  containing  RFI , however,  the  m/a  ratio  sharply  in- 
creases. 

The  BRISA  algorithm  calculates  the  m/o  ratio  for  each  Doppler. 

For  those  Dopplers  where  m/a  is  below  a threshold  (say,  5)  the  algorithm 
presumes  that  only  noise  or  clutter  is  present,  and  the  standard  ERS 
normalization  is  performed.  For  Dopplers  with  m/a  greater  than  5,  BRISA 
first  subtracts  the  mean  amplitude  at  that  Doppler  (m)  from  all  amplitudes 
at  that  Doppler.  Since  this  produces  both  positive  and  negative  values 
(a  strong  target  might  have  a large  negative  value,  depending  upon  its 
phase  relationship  to  the  RFI)  , the  absolute  value  of  each  amplitude  is 
taken  after  subtraction  of  the  mean.  Finally,  the  standard  ERS  normaliza- 
tion is  then  performed.  The  results  of  this  process  are  shown  in  Fig- 
ure 4(c). 

3.  A Test  of  BRISA 

For  the  purpose  of  this  brief  description  it  is  appropriate  to 
provide  an  example  of  the  performance  of  BRISA.  Radar  data  were  collected 
using  the  WARF  radar.  During  data  recording,  a local  RFI  source  (with  ap- 
proximately 25  dB  SNR)  was  switched  on  and  its  frequency  was  controlled  to 
place  the  RFI  on  the  radar  output  at  the  same  Doppler  as  an  observed  target 
of  opportunity.  These  data  were  recorded  to  be  used  in  later  analysis  of 
different  normalization  techniques.  Figure  6 displays  a 15-minute  data 
sequence  in  the  range-versus-time  display  (nested  Doppler)  format  to  be 
used  at  the  ERS.  Standard  ERS  normalization  has  been  performed. 

In  Figure  7 the  same  data  have  been  reprocessed  using  BRISA. 

The  increase  in  track  detectability  is  apparent.  In  Figure  8 the  target's 
SNR  ver  !us  time  is  plotted  for  the  standard  ERS  normalization  (Figure  6) 
and  for  BRISA  normalization  (Figure  7).  Figure  9 replots  the  data  of 
Figure  8 in  histogram  fashion.  These  results  demonstrate  that  application 
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, . . . j Dcn  PRnrF'iSED  USING  ERS  NORMALIZATION 

RADAR  DATA  (with  targets  and  RFII  PHOLtbbtu 


FIGURE  7 RADAR  DATA  (with  targets  and  RFI)  PROCESSED  USING  BRISA 


of  BRISA  can  make  significant  difference  for  target  detectability  in  the 
presence  of  RFI . 


4.  Discussion 

Several  observations  made  during  the  development  of  BRISA  are 
listed  below: 

• The  amount  of  RFI  that  can  be  suppressed  by  BRISA  is 
limited  by  the  accuracy  with  which  the  receiver  pass- 
band  is  known.  (WARF  results  indicate  that  30-dB 
suppression  is  easily  achievable.) 

• Degradation  of  the  algorithm  performance  was  observed 
as  the  RFI  source  was  positioned  toward  the  edges  of 
swept  bandwidth.  (The  controlled  RFI  source  used  in 
the  data  sequence  shown  above  was  placed  approximately 
at  the  middle  of  the  swept  bandwidth.  Over  20  dB  sup- 
pression was  achieved  for  more  than  504  of  the  swept 
bandwidth . ) 

• BRISA  may  need  modification  for  use  with  radar  wave- 
form parameters  that  produce  a high  RFI  time-bandwidth 
product  (radar  processor  IF  bandwidth  times  the  length 
of  time  the  RFI  is  in  that  bandwidth  each  w.'iveform 
sweep).  The  parameters  used  in  collecting  the  data  of 
Figures  6 and  7 have  an  RFI  time-bandwidth  product  of 
1.  However,  no  noticeable  degradation  in  suppression 
performance  was  observed  when  a time-bandwidth  product 
of  10  was  used. 

• BRISA  worked  as  well  on  real  RFI  sources  (i.e. , sky- 
wave  signals)  as  it  did  on  the  local  RFI  source  used 
in  the  above  measurements. 

• BRISA  can  potentially  work  on  all  components  of  RFI. 

BRISA  performance  is  not  limited  by  a spread-Doppler 
RFI  source  as  long  as  its  spectrum  remains  constant 
as  the  RFI  is  swept  through  the  receiver  passband. 

In  summary,  BRISA  promises  to  be  an  attractive  solution  to  the 
RFI  problem  foreseen  for  the  ERS.  As  has  been  pointed  out  by  MITRE,  the 
typical  RFI  environment  will  consist  of  several  (3  to  5)  carrier  sources, 
each  having  approximately  a 20-dB  SNR.  With  the  current  normalization, 
such  RFI  will  result  in  SNR  degradations  of  20  dB  for  several  target 
speeds.  Application  of  BRISA  could  reduce  this  degradation  to  less  than 
3 dB  for  each  of  these  speeds.  Coupled  with  the  frequency-stepping  tech- 
nique proposed  by  MITRE,  application  of  BRISA  should  significantly  improve 
ERS  performance  during  the  difficult  nighttime  RFI  environment. 
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c. 


Noise  Model  and  Display  Evaluation 


In  1976  SRI  developed  an  analytical  model  that  simulates  OTH-radar 
target  and  noise  data.®  During  the  current  reporting  period  this  data 
model  was  implemented  on  the  Hewlett-Packard  2100  minicomputer  facility 
located  at  the  Remote  Measurements  Laboratory  of  SRI.  This  facility  has 
the  capability  of  generating  both  the  WARE  animated  display  and  the 
ERS  gray-scale  display.  Thus,  computer  implementation  of  the  data  model 
provides  a controlled  environment  for  the  test,  evaluation,  and  comparison 
of  these  display  formats. 


1.  A Teat  of  Model  Realism 

The  analytical  data  model  was  developed  for  the  purpose  of  pro- 
viding realistic  targets  and  noise.  To  accomplish  this,  the  model  was 
based  on  target  and  noise  statistics  extracted  from  real  OTH-radar  data 
collected  at  WARE.  To  determine  how  well  the  model  simulates  actual  OTH 
radar  data,  a test  run  was  conducted  to  recreate  a specific  target  track 
obtained  at  WARE.  Figure  10  presents  a range  (with  nested  Doppler) -vs- 
timc  display  of  WARE  radar  data  showing  two  aircraft  tracks.  Figure  11 
presents  a recreation  of  Figure  10  using  the  OTH  data  model.  The  Dopplers 
of  the  recreated  tracks  do  not  precisely  coincide  with  those  of  the  orig- 
inal tracks  because,  for  simplicity,  the  model  assumes  an  operating  fre- 
quency of  15  MHz  while  the  data  were  actually  collected  at  a higher 
frequency.  Other  noticeable  differences  include: 

• Target  f ading--Dwell-to-dwell  SNR  fluctuations  are  part 
of  the  data  model.  However,  long-term  fading  (ove. 
several  dwells)  such  as  can  be  seen  in  Figure  IJ  is  not 
included  in  the  model. 

• Target  signature  changes — A target  signature  function 
is  included  in  the  model  to  simulate  range  and/or  Dop- 
pler spreading  of  targets.  This  signature  function  is 
restricted  to  remain  constant  during  the  progress  of 
each  track.  However,  in  the  case  of  actual  targets, 
signatures  can  change  as  the  target  proceeds  (e.g. , 
range  spreading  decreases  as  targets  increase  in  range 
in  Figure  10) . 

• Meteors--The  model  has  the  capability  of  including  ran- 
domly occurring  meteors  at  an  expected  average  rate  of 
occurrence.  (Meteors  appear  s the  small  dark  rectangles 
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SLANT-RANGE  OFFSET 


on  the  displays.)  However,  a comparison  of  Figures  10 
and  11  indicates  that  the  signature  function  and  ampli- 
tude distribution  function  used  to  simulate  meteors 
could  be  improved. 

For  the  most  part  these  differences  are  minor.  The  major  char- 
acteristics of  OTH  targets  and  noise  are  reasonably  well  reproduced  using 
the  data  model.  Thus  we  feel  justified  in  using  this  model  for  compara- 
tive evaluation  of  detection  and  tracking  systems.  The  first  use  of  this 
model  was  directed  toward  determining  the  detection  sensitivity  of  the 
ERS  display. 

2.  Display  Sensitivity 

The  data  model  is  able  to  create  displays  of  realistic  targets 
with  specified  average  SNR  and  fluctuation  levels,  and  is  thus  ideally 
suited  for  testing  display  sensitivity.  For  an  initial  test,  a display 
of  range  (nested  Doppler)  vs  time  containing  targets  of  0 to  20  dB  SNR  wa 
generated.  Figures  12(a),  13(a),  and  14(a)  reproduce  the  displays,  while 
Figures  12(b),  13(b),  and  14(b)  indicate  track  positions,  SNRs,  and  fluc- 
tuation rates. 

One  useful  measure  of  display  sensitivity  can  be  defined  as  the 
SNR  required  to  obtain  a 507.  detection  rate.  From  analysis  of  these  data 
it  appears  that  the  display  sensitivity  lies  between  6 and  9 dB. 

The  test  described  above  served  to  determine  the  range  of  SNRs 
that  have  marginal  detection  rates.  In  the  coming  months,  additional 
simulations  will  be  generated  to  obtain  a more  precise  estimate  of  dis- 
play sensitivity.  These  simulations  will  also  be  used  to  evaluate  the 
display  of  Doppler  (nested  range)  vs  time  as  well  as  the  WARF  animated 
display.  Finally,  the  data  model  will  be  used  to  test  the  automatic 
detection  and  tracking  algorithm  proposed  for  the  ERS.  Besides  deter- 
mining the  sensitivity  of  this  algorithm,  it  is  hoped  that  we  can  gain 
insight  into  appropriate  methods  of  interaction  between  the  manual  and 
automatic  detection/tracking  process. 
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D. 


Generalized  Sldelobe  Canceller 


The  University  of  Colorado  (UC)  under  subcontract  to  SRI  has  devel- 
oped a new  adaptive  noise-cancelling  scheme  referred  to  as  a Generalized 
Sidelobe  Canceller  (GSC).  This  approach  to  adaptive  signal  processing 
may  have  application  to  OTH  radars  operating  in  the  presence  of  auroral 
clutter.  This  brief  description  of  the  GSC  processor  will  focus  on  its 
conceptual  attributes;  an  in-depth  analysis  will  be  presented  in  a future 
UC  report.  In  the  meantime,  while  the  GSC  structure  is  under  development 
at  Colorado,  a version  of  this  processor  has  been  implemented  in  software 
at  SRI  for  off-line  processing  of  WARF  data  and  for  preliminary  performance 

evaluation. 

Figure  15  is  a simplified  diagram  of  the  GSC  processor.  A multidi- 
mensional input  signal  X as  a function  of  time  K in  the  present  applica- 
tion represents  the  signals  '■eceived  at  an  array  of  spaced  antennas.  This 
array  or  vector  of  input  si  Is  is  fed  through  two  signal  paths  simul- 
taneously, referred  to  as  the  upper  and  lower  paths  of  Figure  15.  The 
function  of  the  lower  path  is  first  to  apply  a spatial  preprocessing 
filter  W to  the  data  in  order  to  remove  the  desired  signal  from  the  data. 
The  residue  at  the  output  is  termed  noise  or  unwanted  energy.  These  data 
are  fed  through  an  adaptive  feedback  loop  , which  in  general  involves 
both  spatial  and  time-domain  filtering,  to  generate  a beamformed  output 

v,(k). 

The  processing  functions  of  the  upper  loop  are  a beamforming  opera- 
W followed  by  a feedback  tapped-delay-llne  W^.  The  output  of  the  upper 
path  Y^(k)  is  subtracted  from  the  lower  path  output  to  generate  the 

final  GSC  output  signal,  Y^(k) . Note  that  Y^{k)  serves  as  the  feedback 
error  signal  to  update  the  adaptive  weights  in  both  the  upper  and  lower 
paths.  The  weight-updating  algorithm  for  the  lower  path  is  the  least- 
mean-squares  (LMS)  algorithm,""  while  weight  updating  for  the  upper  path 
may  take  on  several  forms. 

The  lower  path  seir^es  as  an  error  beam  (or  perturbation  signal)  so 
that  if  unwanted  signals  are  present,  the  process  of  subtraction  Y^(k) 
from  Y '(k)  will  remove  such  signals  from  the  final  output.  If  no  such 
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FIGURE  15  GENERALIZED  SIDELOBE  CANCELLER 


signals  are  present,  the  error  beam  automatically  shuts  down  and  the  uppe 
path  continues  to  process  the  input  signals  in  a manner  similar  to  that 
of  a conventional  beamformer. 

One  of  the  most  promising  aspects  of  this  processor  structure  is 
that  it  permits  the  realization  of  a fairly  general  class  of  adaptive 
beamforming  algorithms  with  a single  unifying  analytical  approach.  This 
feature  greatly  enhances  the  ability  to  compare  the  effect  of  different 
approaches  to  adaptive  beamforming,  something  that  has  not  been  readily 
possible  in  previous  work.  For  example,  it  is  now  quite  easy  to  make 
this  beamformtng  structure  realize  the  Frost  algorithm  or  the  Griffiths 
algorithm,  two  previously  studied  approaches  to  adaptive  beamforming. 
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Details  of  these  versions  of  the  GSC  will  be  presented  in  a technical 
report  devoted  to  this  subject. 

One  current  area  of  interest  at  UC  is  the  form  that  the  spatial  pre- 
processing matrix  Hg  should  take  in  order  to  most  effectively  remove  the 
want'-d  signal  from  the  error  beam  One  of  the  simplest  versions  is  to  use 
this  filter  to  combine  the  array  outputs  in  pairwise  monopulse  fashion-- 
that  is,  element  1 is  subtracted  from  element  2,  element  2 is  subtracted 
'rom  element  3,  and  so  on.  A plane-wave  signal  in  phase  across  ill  array 
elements  will  be  removed  by  this  process,  leaving  the  so-called  unwanted 
signals  or  noise.  At  UC  Griffiths  has  shown  how  multielement  generaliza- 
tions of  this  monopulse- like  spatial  preprocessor  may  be  realized  using 
Wals^  functions. 


The  GSC  has  been  implemented  in  minicomputer  software  at  SRI  for  off- 
line processing  of  WARF  data.  This  implementation  is  quite  general  and 
flexible,  to  facilitate  comparison  of  different  realizations  of  the  various 
component  parts  of  the  GSC.  The  present  configuration  has  the  following 
options: 


• weights 

• ^s 


Uniform  or  Lilph 

Structured  to  realize  tither  Frost 
or  Griffiths  algorithms  with  8 
elements  and  5 taps 

Pairwise  monopulse,  Walsh  functions, 
or  manually  selected  at  program 
initiation 


• Elemental  power  Signals  at  output  of  Hg  have 

equalization  integrate-and-dump  normalization 

with  operator-selectable  exponen- 
tial time  constant  (usually  100 
samples) 

• Range-Doppler  analysis  YoC’*^)  is  applied  to  a 2-s  (4096- 

sample)  range-Doppler 

• Output  oisplay  Either  range-Do.»iAer  map  or  Doppler 

power  spectrum. 


Debugging  of  this  processor  was  achieved  in  May  through  the  use  of 
taped  artificial  data  designed  to  simulate  both  desired  and  unwanted  sig- 
nals arriving  at  an  8-element  array  with  320-m  spacing.  Based  or  nese 
successful  tests,  the  GSC  was  used  in  a preliminary  test  of  beamfo^mer 
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performance  against  azimuthally  spread  ground  clutter;  the  results  of  this 
test  are  described  in  Section  II-E. 

Two  tests  of  the  GSC  processor  were  carried  out  with  the  aid  of  simu- 
^-ited  signals.  In  the  first  test,  two  sinusoidal  signals  were  injected 
into  the  eight  input  channels,  one  signal  in-phase  across  the  eight  chan- 
nels (in-beam) , and  the  other  signal  with  .a  linear  phase  shift  to  make 
the  signal  appear  l.S”  away  from  the  beam  pointing  direction.  Gaussian 
noise  was  added  so  tnat  the  bNR  of  the  in-beam  signal  was  30  dE  when  mea- 
sured after  beamforming  and  a 4096-sample  coherent  integration  (=“2  s). 

The  out-of-beam  signal  would  have  had  a 50-dB  SNR  if  it  were  in  the  beam. 
The  desired  (in-beam)  signal  was  at  a frequency  of  420  Hz,  and  the  un- 
wanted signal  at  400  Hz.  In  the  second  test  all  conditions  were  identical 
except  that  both  signals  were  given  element-to-element  amplitude  and 
phase  randomization  (107«  and  15**  respectively).  This  randomization  was 
different  for  each  of  the  two  signals  and  it  was  fixed  in  time  throughout 
the  4096  samples.  Its  purpose  was  to  simulate  typical  departures  f.om 
equal  amplitude  and  linear  phase  front  that  are  seen  in  real  signals. 

These  signals  were  applied  to  the  GSC  using  both  the  Griffiths  and 
Frost  beamforming  algorithms  with  variable  gain  (a)  in  the  weight- 

recursion  loops.  Conventional  weights  W were  25-dB  Dolph  coefficients; 

“"C 

the  W preprocessing  filter  was  the  pairwise-monopulse  matrix--!. e. , 
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The  time  constant  of  the  individual  channel  gain  normalization  was  100 
samples. 
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Figure  16  presents  the  Frost  algorithm  results.  SNR  is  plotted  for 
loop  gain  o for  both  desired  and  unwanted  signals  and  for  both  ideal  and 
randomized  amplitudes  and  phases.  Desired  signal  SNR  remains  near  30  dB 
at  small  values  of  a but  begins  to  decrease  sharply  for  a > 0.1.  This 
reduction  in  SNR  is  expected  and  is  caused  by  adaptation  noise.  Cancel- 
lation of  the  unwanted  signal  is  complete  at  o?  =•  0.001  for  ideal  signals 
and  a =“  0.005  for  the  randomized  signal.  These  results  indicate  that 
the  correct  gain  has  a fairly  wide  range  of  acceptable  values  centered 
around  0.01. 

Corresponding  plots  for  processing  via  the  Griffiths  algorithm  are 
given  in  Figure  17.  The  overall  trends  are  similar  to  those  for  Frost 
processing  except  that  a slightly  larger  value  of  o?  is  required  for  the 
Griffiths  processing  to  achieve  the  same  amount  of  rejection.  Once 
again,  higher  a is  required  to  remove  the  unwanted  signal  for  the  ran- 
domized signal  case.  One  notable  feature  of  Griffiths  beamforming  is 
that  SNRs  for  both  signals  appear  larger  than  those  measured  for  the 
conventional  beamformer.  This  behavior  is  caused  by  the  tapped-delay- 
line  filtering,  which  is  an  integral  part  of  the  Griffiths  algorithm 
and  which  produces  filtering  (i.e. , noise  cancellation)  in  the  frequency 
domain.  The  noise- measurement  algorithm  used  for  these  data  computes  an 
rms  value  over  the  whole  frequency  band,  and  the  computed  noise  values 
may  be  decreased  since  some  noise  is  cancelled  during  adaptation.  The 
bias  in  noise  values  could  be  removed  by  measuring  noise  at  the  mid- 
frequencies of  the  receiver  passband  (~i20  Hz) , but  since  both  wanted 
and  unwanted  signals  are  referenced  to  the  same  noise  measurement  it  is 
of  no  consequence. 

These  results  are  in  general  accordance  with  the  expected  performance 
of  the  GSC  for  the  simulated  signals  used  In  the  test,  and  therefore  it 
is  concluded  that  the  SRX  implementation  is  functioning  properly.  Fur- 
thermore, these  results  establish  the  desired  amount  of  feedback  in  the 
weight-recursion  algorithms. 
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SNR — dB  SNR 


LOOP  GAIN  — 0. 

(a)  IDEAL  SIGNALS  (IN-BEAM  SNR  - 30  dB; 
OUT-OF-BEAM  SNR  - 17  dB) 


(IN-BEAM  SNR  - 31  dB;  OUT-OF-BEAM  SNR  - 33  dB) 

FIGURE  16  PERFORMANCE  OF  GSC  PROCESSOR  (Frost  Algorithm)  AGAINST  SIMULATED 
SIGNALS 
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(a)  IDEAL  SIGNALS  (IN-BEAM  SNR  - 30  dB; 
OUT-OE-BEAM  SNR  - 17  dB) 


LOOP  GAIN  — a 

(b)  AMPLITUDE  AND  PHASE  RANDOMIZED  SIGNALS 
(IN-BEAM  SNR  - 31  dB;  OUT-OF-BEAM  SNR  - 33  dB) 

F(GURE  17  PERFORMANCE  OF  GSC  PROCESSOR  (Griffitlis  A(gorit)im)  AGAINST  SIMULATED 
SIGNALS 


cl.  Spread  Clutter  Effects  on  Adaptive  Beamformers 

One  of  the  primary  reasons  for  applying  adaptive  beamforming  tech- 
niques to  HF  radar  is  to  exploit  their  potential  for  rejecting  spread- 
azimuth  and  spread-Doppler  clutte,;  such  as  that  produced  by  scattering 
from  auroral  irregularities.  While  it  is  relatively  easy  to  implement 
point-source  interferers  for  performance  testing  of  a candidate  beam- 
forming technique,  it  is  not  convenient  to  generate  a realistic  spread- 
azimuth  interference  source.  Computer  simulation  studies  have  been  under- 
taken at  SRI  to  analyze  the  degree  to  which  various  methods  of  adaptive 
beamforming  could  be  expected  to  work  against  such  unwanted  energy.  Yet 
the  simulation  results  are  not  fully  convincing  in  their  ability  to  ac- 
count for  all  the  aspects  of  the  actual  HF  environment  and  the  achievable 
performance  of  HF  equipment.  Experimental  results  are  highly  desirable. 

Consideration  has  been  given  at  SRI  to  determining  how  the  WARF 
adaptive  beamforming  capabilities  could  be  used  to  test  ways  to  reject 
spread-azimuth  interference.  Possible  vjays  to  address  this  problem  ex- 
perimentally include  the  following: 

(1)  Use  a separate  transmitter--The  WARF  could  be  configured 
to  look  at  desired  targets  such  as  a repeater  while  a 
remote  transmitter  operating  on  the  same  frequency  is 
located  to  make  energy  arrive  through  the  sidelobes  of 
the  receive  beam.  If  this  out-of-beam  signal  propagates 
via  a 2-hop  path,  then  some  azimuthal  and  Doppler  spread- 
ing of  the  unwanted  signal  will  result.  Depending  pri- 
marily on  the  transmit  beamwidth  and  the  terrain  roughness 
at  the  midpath  point,  the  azimuthal  spreading  of  the  re- 
mote transmitter  signal  could  amount  to  several  degrees.^ 

(2)  Mis- steer  the  WARF  transmit  beam--The  log  periodic  array 
at  Lost  Hills  has  a 6°  azimuthal  beamwidth  (at  15  MHz)  with 
average  sidelobes  20  to  25  dB  down.  By  the  simple  expe- 
dient of  shifting  the  transmit  steer  direction  away  from 
the  direction  to  the  desired  signal , a clutter  distribution 
with  significant  out-of-beam  components  can  be  produced  at 
the  Los  Banos  receive  site.  The  out-of-beam  clutter  has 
about  6°  azimuthal  width  and  a Doppler  spectrum  like  that 
of  normal  land  or  sea  backscatter. 

(3)  Gather  data  in  the  presence  of  spread  ionization--Even  at 
the  moderate  latitude  of  WARF  (55°  magnetic),  ionospheric 
phenomena  similar  to  those  observed  in  the  auroral  regions 
are  present.  Scatter  from  spread-F  ionization  is  observed 
frequently;  and  reflections  from  these  irregularities 
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produce  spreading  in  Doppler  and  in  azimuth.  Data  col- 
lection in  the  presence  of  such  phenomena  could  emulate 
almost  exactly  the  eventual  operational  environment  of 
interest  to  the  Air  Force, 

(4)  Use  the  Boulder  ionospheric  heater--IonosDheric  heaters 
generate  field-aligned  ionospheric  irregularities  very 
much  like  those  found  in  auroral  regions.  A critical 
and  unique  aspect  of  these  artifically  induced  irregu- 
larities (and  therefore  sources  of  clutter)  is  that  they 
are  controllable--i.e. , they  may  be  turned  on  and  off  at 
will.  This  property  makes  heater-induced  spread  clutter 
more  desirable  than  naturally  occurring  spread  clutter  at 
least  for  purposes  of  evaluating  beamforming  techniques. 

At  the  present  time  only  the  first  alternative — a separate  transmitter- 
has  been  rejected,  principally  because  of  the  anticipated  difficulties  in 
frequency  management  for  simultaneous  one-hop  and  two-hop  paths.  Mis- 
steering  of  the  transmit  beam,  while  not  fully  realistic,  is  inexpensive 
and  simple  to  test.  Some  WARF  data  collected  with  this  configuration  are 
described  in  Section  II -F,  below.  The  third  and  fourth  options  of  collect- 
ing data  in  the  presence  of  natural  or  artificial  spread  clutter  may  be 
tested  in  FY78. 

F.  Mis-Steered  Transmitter  Beam 

Steering  the  transmit  beam  away  from  the  desired  target  (the  Los 
Lunas  New  Mexico  repeater  in  this  case)  produces  substantial  ground 
clutter,  which  is  received  through  sidelobe  angles.  This  clutter  has 
azimuthal  extent  at  least  as  large  as  the  transmit  beamwidth  (6°),  and 
in  the  frequency  domain  consists  of  about  1-Hz-wide  components  spaced  at 
the  waveform  repetition  rate.  This  situation  provides  a preliminary  test 
of  beamformer  capability  to  reject  spread  clutter. 

The  WARF  was  operated  as  in  earlier  adaptive  beamforming  experi- 
ments^— that  is,  essentially  in  an  aircraft  detection  mode  with  the  ex- 
ception of  the  eight-channel  coherent  receiving  system.  In  Phase  1 of 
the  test  the  transmitted  beam  was  pointed  at  the  repeater  as  accurately 
as  possible  while  several  data  fries  of  the  receiver  outputs  were  re- 
corded digitally.  Sequential  files  were  recorded  as  the  conventional 
array  pointing  direction  was  stepped  ±1/4°  from  the  nominal  azimuth  to 
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the  New  Mexico  repeater;  this  maneuver  is  included  to  ensure  that  the 
repeater  echo  is  closely  aligned  to  the  pointing  direction  in  at  least 
one  of  three  sequential  data  files.  Phase  2 of  the  test  repeated  the 
sequential  stepping  of  the  receive  look  direction  while  the  transmit 
steer  was  misaligned  hv  12’.  The  sketch  in  Figure  18  illustrates  the 
the  time  sequence  of  the  test. 


RECEIVER 


SEQUENTIAL 
STEPPING 
OF  RECEIVE 
LOOK  DIRECTION 


FIGURE  18  GEOMETRY  OF  TRANSMITTER  MIS-STEERING  TEST  (not  to  scale). 

At  15  MHz,  receive  beamwidth  is  0.5°  and  transmit  beamwidth  is  6 . 


Over  a period  of  6 minutes,  72  two-second  coherent-integration  pe- 
riods were  acquired,  36  with  the  transmitter  beam  pointing  as  close  to 
the  repeater  as  possible  and  36  with  the  beam  misaligned  by  12  . The 
repeater  gain  was  increased  by  30  dB  during  the  misaligned  portion  of 
the  test  in  order  to  maintain  repeater  SNR  at  a mean  value  of  25  dB  over 
both  portions  of  the  test.  Clutter-r.o-noise  ratio  (or  subclutter  visi- 
bility) decreased  by  more  than  20  dB  when  the  transmitter  beam  was  mis- 
steered.  Figure  19  presents  typical  samples  of  total  power  spectra  for 
the  two  phases  of  the  test;  these  spectra  were  obtained  using  conventional 
beamforming  and  a 2-s  (or  128-sweep)  coherent  integration  time. 

A version  of  the  Generalized  Sidelobe  Canceller  (GSC)  beamformer  was 
applied  to  these  data.  A matrix  consisting  of  pairwise  monopulse 
weights  (described  in  Section  II-D)  was  used  to  cancel  the  desired  signal 
in  the  error  (perturbation)  beam  processing  patn.  The  Frost  algorithm 
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using  35  weights  (seven  elements,  five  taps)  with  a = 0.1  was  applied  to 
this  error  signal  to  provide  the  error  beamformed  output.  This  output 
was  subtracted  from  the  conventional  beam  and  the  result  underwent  range- 
Doppler  analysis  as  in  the  conventional  case. 

Table  1 summarizes  the  data  analysis  for  the  two  sets  of  data  by 
presenting  average  values  of  SNR  and  signal-to-clutter  ratio  (SCR).  In 
each  case,  signal  is  defined  as  the  peak  signal  amplitude;  clutter  is 
defined  as  the  maximum  value  (over  Doppler)  of  the  rms  range  average  of 
clutter  magnitude  (rms  over  216  km),  and  noise  is  defined  as  the  rms  of 
noise  over  regions  of  range  and  Doppler  sufficiently  removed  from  clutter 
and  targets.  These  quantities  are  measured  after  range-Doppler  analysis, 
and  in  order  to  refer  them  to  the  point  of  beamformer  output  we  subtract 
33  dB  from  SNR  and  20  dB  from  SCR.  Since  processing  gain  affects  both 
conventional  and  adapted  channels  equally,  comparisons  are  made  after 
range-Doppler  analysis.  These  quantities  are  in  fact  measured  from  data 
plots  such  as  those  in  Figure  19, 


Table  1 

DATA  SUMMARY  FOR  MIS -STEERED  TRANSMITTER  BEAM  TEST 


Transmit  Steer 
90° 

Transmit  Steer 
102" 

Conventional  SCR  (dB) 

-32 

-14 

SCR  Improvement  (dB) 

8 

11 

Conventional  SNR  (dB) 

25 

26 

SNR  Improvement  (dB) 

-3.9 

0.3 

In  Table  1,  improvement  in  SNR  and  SCR  refers  to  increases  in  ..nese 
quantities  through  adaptive  beamforming  over  those  achieved  by  conventional 
beamforming.  On  the  average,  the  SCR  improvement  is  8 dB  with  beam  align- 
ment, increasing  to  11  dB  with  beam  misalignment.  The  SCR  improvement 
viith  beam  misalignment  would  in  all  likelihood  be  larger  if  the  conven- 
tional SCR  had  been  constant,  since  it  is  known  that  the  potential  for 
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adaptive  beamformer  improvement  is  larger  for  small  ratios  of  desired 
signal  to  unwanted  signal  at  the  input  to  the  processor.  Thus,  if  con- 
ventional SCR  could  have  been  maintained  at  -32  dB  when  the  beam  was 
mis-steered,  then  SCR  improvement  should  have  increased  more  than  from 
8 to  11  dB,  The  point  here  is  that  even  with  the  presence  of  quasi- 
grating lobes  the  adaptive  beamformer  achieves  cancellation  of  the  spa- 
tially extended  ground  clutter, 

A related  quantity  of  interest  is  SNR,  which  is  expected  to  remain 
constant  if  the  noise  is  spatially  Isotropic,  However,  if  the  repeater 
signal  is  not  precisely  within  the  received  beam,  the  beamformer  may 
partially  reject  the  signal,  and,  furthermore,  with  high  gain  (Of)  in  the 
adaptive  loop  some  adaptation  noise  may  be  added,  SNR  improvement  with 
beam  misalignment  is  virtually  0 dB  as  expected;  however,  with  beam 
alignment  it  is  -4  dB,  Since  the  data  conditions  were  identical  for 
both  data  sets,  it  is  conjectured  that  the  larger  value  of  clutter  (rela- 
tive to  the  signal)  for  the  data  with  transmit  beam  alignment  has  produced 
more  adaptation  noise.  This  source  of  noise  is  redu^ed  with  smaller 
values  of  a;  however,  lower  or  may  possibly  reduce  the  achieved  amount  of 
SCR  improvement. 

This  preliminary  test  of  the  CSC  performance  against  spatially  spread 
clutter  is  not  conclusive,  but  it  is  encouraging.  First,  some  portion  of 
the  clutter  is  either  within  the  beam-pointing  direction  or  within  a 
quasi-grating-lobe  pointing  direction,  and  therefore  an  adaptive  beamformer 
can  only  partially  cancel  the  clutter.  Thus,  this  test  cannot  be  ex- 
pected to  measure  the  achievable  amount  of  cancellation  of  out-of-beam 
spatially  spread  clutter.  Second,  the  frequency  spread  of  this  type  of 
clutter  consists  of  discrete  frequencies  spaced  at  the  waveform  repetition 
rate  with  minor  Doppler  spreading  (1  to  2 Hz)  about  these  frequencies. 
Ultimately  we  would  like  to  measure  performance  against  clutter  with 
larger  Doppler  frequency  spreading.  On  the  other  hand,  this  result  does 
show  that  appreciable  cancellation  of  azimuthally  spread  clutter  is 
achievable.  Future  experimentation  will  seek  to  determine  expected  per- 
formance levels  in  a more  realistic  experimental  environment. 
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Ill  CONSULTING  ACTIVITIES 


During  the  six  months  covered  by  this  report  SRI  provided  a variety 
of  services  requested  by  the  SPO  under  the  consulting  portion  of  this 
contract.  The  following  is  a brief  description  of  each  of  the  major 
consulting  tasks: 

(1)  Cost-Savings  Ef fort--Because  of  growing  costs,  the  Prototype 
Radar  System  (PRS)  design  was  substantially  modified  to  the 
current  Experimental  Radar  System  (ERS)  design.  SRI  Inter- 
national supported  this  effort  by  participating  in  meetings  at 
the  SPO,  providing  timely  data  to  the  SPO  for  high-level  brief- 
ings, and  presenting  technical  opinions  in  response  to  requests 
from  the  Office  of  the  Assistant  Secretary  of  the  Air  Force  for 
Research  and  Development. 

(2)  Long-Path  Soundings--At  the  request  of  ESD,  SRI  provided  ap- 
proximately 200  hours  of  sounding  transmissions  from  the  WARF 
Lost  Hills  transmit  facility.  These  soundings,  requested  by 
Dr.  Elkins  of  RADC,  were  received  on  RADC  equipment  at  Ava , 

New  York  and  later  overseas. 

(3)  Transponder  Design — It  had  been  proposed  that  SRI  supply  the 
Air  Force  with  transponders  for  use  in  the  ERS  coverage.  Pre- 
liminary design  and  cost  estimates  were  developed  during  this 
reporting  period. 

(4)  Environmental  Assessment  System--SRI  developed  a technical  con- 
cept for  augmenting  the  ERS  with  an  independently  operated  en- 
vironmental assessment  system  (EAS) . This  would  be  a 
minicomputer-based  system  designed  to  improve  ionospheric 
propagation  management  capabilities  at  the  ERS. 
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Appendix  A 

ARCHITECTURE  FOR  A HARDWARE  BEAMFORMER 


As  a part  of  the  ongoing  analysis  of  adaptive  beamformers  at  SRI,  a 
study  has  been  initiated  to  define  how  an  on-line  beatnformer  might 
e ’entually  be  best  implemented.  During  this  reporting  period  the  struc- 
ture of  the  CSC  was  used  as  a candidate  for  possible  on-line  implementa- 
tion in  order  to  obtain  a measure  of  cost,  speed,  and  performance  of 
such  a device  when  integrated  into  the  operation  of  an  OTH  radar. 

The  particular  structure  of  the  generalized  sidelobe  cancelle*'  has 
led  to  the  suggested  design  of  a multipler-oriented  processor  shown  in 
Figure  A-1.  This  design  takes  advantage  of  the  fact  that  many  of  the 
algorithm  operations  are  of  the  form  of  a weight  term  times  a data  or 
output  term.  This  simplifies  the  splitting  of  memory  into  two  parts, 
one  a data  memory  and  one  a weight  memory,  which  decreases  the  access 
time  for  the  multiple!  and  multiplicand.  Also  some  of  the  equations  in- 
volve accumulating  the  products  of  data  and  weights  for  up  to  35  opera- 
tions, with  the  results  stored  in  either  memory.  This  fact  led  to  the 
design  of  the  product  summation  unit  as  an  adder-subtractor-accumulator 
with  access  to  either  memory. 

Some  of  the  algorithms  involve  summing  data  in  memory.  This  could 
be  accomplished  by  multiplying  by  1 if  integer  arithmetic  is  assumed,  or 
by  the  largest  positive  fraction  if  fractional  arithmetic  is  a.ssumed,  -^r 
by  -1  and  subtracting  the  result  if  integer  arithmetic  is  assumed.  An 
easier  method  is  to  allow  the  memories  to  have  direct  access  to  the  ac- 
cumulator through  a bus  and  selection  arrangement.  This  method  permits 
easier  and  quicker  summation  of  nonmultiplied  terms. 

In  order  to  prevent  any  loss  of  accuracy,  the  accumulator  must  be 
designed  to  handle  the  summation  of  up  to  35  products,  each  product  hav- 
ing 23-bit  resolution,  resulting  in  29  bits  of  storage  for  the  accumula- 
tor. The  output  to  memory  can  contain  only  12  bits,  since  integer 
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FIGURE  A-1  ARCHITECTURE  FOR  AN  OPERATIONAL  GENERALIZED  SIDELOBE  CANCELLER 


arithmetic  is  assumed,  and  must  therefore  be  truncated  or  rounded.  This 
illustrates  the  usefulness  and  necessity  of  the  floating-point  approach 
in  minimizing  the  errors  due  to  truncation,  rounding,  or  an  arbitrary 
selection  of  bits. 

The  control  of  the  multiplier,  accumulator,  and  data  buses  would  be 
done  with  Read  Only  Memory  (ROM).  The  ROM,  addressed  by  a counter, 
would  determine  the  addresses  for  each  memory  and  would  control  the  bus 
access  of  each  memory  or  multiplier  output.  Addresses  for  data  storage 
as  well  as  read-write  controls  are  also  generated.  Control  of  the  accumu- 
lator operation  by  the  ROM  provides  clear  accumulator,  accumulator-to- 
output-bus,  sum-input,  and  subtract-input  signals. 

The  ROM  would  also  control  the  data  output  bus  and  peripherals  such 
as  the  analog-to-digital  converter  (aDC)  interface,  dirplay  interface, 
and  interface  to  the  main  processor  computer.  The  ADC  is  a stand-alone 
unit  controlled  by  its  own  microprocessor  or  microprogrammed  controller. 
This  peripheral  would  acquire  a set  of  eight  data  samples,  perform  any 
necessary  digital  filtering,  and  signal  the  GSC  controller  that  a set  of 
data  samples  is  ready.  The  computer  interface  communicates  to  the  system 
computer  and  transfers  the  adapted  beamformed  data  to  the  main  computer 
for  further  analysis  (such  as  range-Doppler  processing)  and  output  to  a 
display.  There  would  also  be  a transfer  of  control  information  or  data 
from  the  main  computer  to  the  GSC.  This  would  enable  the  main  computer 
to  alter  or  adjust  the  adaptive  algorithms  or  use  the  GSC  as  an  external 

data  processing  unit. 

By  examining  the  structure  of  the  GSC  and  the  equations  required, 
it  cun  be  seen  that  each  data  set  (eight  inputs  from  the  ADC)  requires 
195  multiplications  and  196  summations  or  subtractions.  Since  this 
arithmetic  must  be  executed  in  less  than  500  |as , a fast  controller  is 
needed  to  address  memory,  control  the  accumulator  and  bus,  and  control 
the  peripherals.  The  ROM  can  do  this  because  of  its  rapid  access  times. 
The  complexity  of  the  operations  and  the  desire  to  modify  the  algorithm 
operation  in  real  time  suggest  that  a microprocessor-based  controller 
should  be  used.  However,  most  MOS  microprocessors  are  too  slow  for  this 
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application,  having  instruction  times  of  1 ns  or  more.  In  contrast,  a 
bipolar  microprocessor  slice  or  microprogram  sequencer  could  be  used  to 
control  the  GSC.  This  alternative  would  permit  adjustment  of  parameters 
to  optimize  GSC  performance  according  to  operational  needs. 
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Appendix  B 

GENERALIZED  SIDELOBE-CANCELLER-SYSTEM  MULTIPLIERS 
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Appendix  B 


GENERALIZED  SIDELOBE-CANCELLER-SYSTEM  MULTIPLIERS 

The  multiplier  used  in  the  adaptive  processor  must  be  able  to  form 
23-bit  products  from  12-bit  two ' s-complement  input  numbers.  The  number 
of  multiplications  required  in  each  data  update  period  (500  |J.s)  is  such 
that  multiplier  speed  must  be  2 p,s  or  less. 

Several  multiplying  schemes  and  implementations  were  studied*  and 
compared  on  the  basis  of  speed,  power  consumption,  total  size  of  multi- 
plier, and  total  costs.  For  comparison,  methods  that  do  not  use  two's- 
complement  arithmetic  or  that  lack  sufficient  precision  are  included. 
Table  B-1  lists  the  costs  of  space  requirements  (board  layout  costs, 
physical  size  of  integrated  circuits),  components,  power,  and  total  size 
for  several  multiplying  implementations,  while  Table  B-2  lists  the  per- 
formance of  these  implementations. 

Studies  are  continuing  to  determine  the  necessary  accuracies  needed 
in  the  GSC  and  to  determine  the  usefulness  of  a floating-point  system. 


* 

"A  Survey  of  Digital  Multiplying  Techniques,"  by  Philip  C.  Evans.  Un- 
published internal  report  for  A.  M.  Peterson  and  T.  W.  Washburn  (June 
1977). 


£Axm  mi  mmi 
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Table  B-1 


CIRCUIT  COSTS 
(Dollars) 


Design 

Requirements 

Space 

Integrated 

Circuit 

Power 

Shaw's  Serial-Parallel  (S-P) 

$23 

$ 53 

$11 

$ 87 

Robertson's  S-P 

24 

38 

9 

71 

Booth's  25LS14 

19 

80 

15 

114 

AMD  2505 

36 

270 

41 

347 

TI  74LS274 

33 

75* 

15 

123* 

TI  74LS261 

41 

106* 

23 

170* 

MM  67558  (4  units) 

22 

418* 

22 

462* 

MM  67558  (1  unit) 

24 

137* 

12 

173* 

TRW  MPY  12AJ 

4 

18 

187 

Squaring  circuit 

28 

113 

9 

150* 

Approximation'^ 

15 

24 

7 

46 

* 

Estimates. 

^Insufficient  precision. 


Table  B-2 

MULTIPLIER  PERFORMANCE 


Design 

Multiplication 

Speed 

(M,s) 

Power 

(W) 

Speed-Power 

(W-M,s) 

S peed -Power- Cost 
($-W-M.s) 

Shaw's  S-P 

1.430 

2.14 

3.06 

266 

Robertson's  S-P 

0.960 

1.68 

1.61 

114 

Booth's  25LS14 

1.000 

2.92 

2.92 

333 

AMD  2505 

0.205 

8.10 

1.66 

576 

TI  74LS274 

0.130 

2.93 

0.38 

47 

TI  74LS261 

0.130 

4.49 

0.58 

99 

MM  67558  (4  units) 

0.140 

4.48 

0.67 

310 

MM  67558  (1  unit) 

0.770 

2.32 

1.79 

310 

TRW  MPY  12AJ 

0.193 

3.50 

0.68 

127 

Squaring  circuit 

0.225 

1.86 

0.42 

63 

Approximation 

0.200 

1.30 

0.26 

12 
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